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Abstract. We study singularity formation in nonlinear differential equations 
of order m ^ 2, j/( m > = A(x _1 , j/). We assume A is analytic at (0, 0) and 
d y A(0,0) = A ^ (say, A = (-l) m ). If m = 1 we assume A(0, ■) is mero- 
morphic and nonlinear. If m = 2, we assume A(0, •) is analytic except for 
isolated singularities, and also that |<E»(s)| 1 / 2 cJ|s| < oo along some path 
avoiding the zeros and singularities of <!>, where Q(s) = f A(0, r)dr. Let 
H a = {z : \z\ > a > 0, arg(^) G (-a, a)}. 

If the Stokes constant S + associated to R + is nonzero, we show that all y 
such that liniaj—^+oo y(x) = are singular at 27ri-quasiperiodic arrays of points 
near iR+. The array location determines and is determined by S + . Such 
settings include the Painleve equations Pj and Pji. If S + = 0, then there 
is exactly one solution yo without singularities in H2 n — e , and yo is entire iff 
y =A(z,0) =0. 

The singularities of y(x) mirror the singularities of the Borel transform of 
its asymptotic expansion, By, a nonlinear analog of Stokes phenomena. If 
m = 1 and A is a nonlinear polynomial with A(z, 0) ^ a similar conclusion 
holds even if A(0, ■) is linear. This follows from the property that if / is super- 
exponentially small along and analytic in H n , then / is superexponcntially 
unbounded in Hn, a consequence of decay estimates of Laplace transforms. 

Compared to [5] this analysis is restricted to first and second order equa- 
tions but shows that singularities always occur, and their type is calculated 
in the polynomial case. Connection to integrability and the Painleve property 
are discussed. 



1. Introduction 

Stokes transitions play an important role in the study of differential equations 
(see pQ, 0], 50j and references therein). In linear differential equations, solutions 
that are small in some sector usually become exponentially large in complementary 
ones. 

There is an interesting analog for analytic nonlinear equations which we study 
for first and second order ones. 

Assume that A{z,y) is analytic at (0,0), and that the following nondegeneracy 
condition holds: d y A(\l : y) := A ^ 0. Let, for m — 1, 2, 

(1) yW=A(l/x,y) 

By a simple change of variable y i— > y + a + bx~ 1 + cx~ 2 we can assume, without loss 
of generality, A(z,0) = o(z 2 ). The conditions in and (which do not require 
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m ^ 2) then apply. By a change of independent variable we can make A = (— l) m . 
Let A[(0) = a. 

There exists a one parameter family of solutions of which decay algebraically 
as x — ► oo ^2] , |3] . All these solutions can be written as generalized Borel summed 
transseries [3] (the notation here is similar to that in 0) 

(2) y(x) = J2 C k e- kx x ak y k (x) = £ C k e~ kx x ak CBy k {x) = £By(x) 

k>0 k>Q 

where y k are formal series in integer powers of x~ x . The functions y k '■= CBy k , 
the generalized Borel sums of y k , are analytic for $l(x) > v > uniformly in k. 
The function series in J5J is absolutely and uniformly convergent for 9?(x) > v 
and defines an analytic function there. Any decaying solution can be represented 
uniquely in terms of Laplace transforms in the upper half plane, y^ — C(By) + 
where + indicates the upper half plane continuation or in terms of the lower half 
plane continuations [3] 

(3) y(x) = C +k e~ kx x ak y+ = C~ h e~ kx x ak y^ 

k>0 k>0 

and we have (0 Eq. (1.27)) 

(4) C+ -C- = S + 

where S + = Si in [3j is the Stokes constant associated to M. + . 

In this paper we prove, under some further assumptions, that all solutions, except 
at most one, that go to zero along some direction form arrays of singularities near 
iR. After normalization of the equation, the array location is given by 

(5) x n = 2mri + a\n(2mri) + InC* - ln£ s + o(l) ±oo) 

for some ^ s ^ 0. (If C _ = (C + = 0), then y has no singularities in the lower 
(upper, respectively) half plane.) In |2j, Eq. (28) we showed that such arrays occur 
provided a solution of an associated equation, JHJ below, is singular. Here we show 
that the solution is singular indeed, under our assumptions. 

There is a correspondence between singularities in inverse Laplace space and 
singularities in the original space, a nonlinear analog of Stokes phenomena: if we 
change variables to make 3^(0,0) = A, then n becomes n/\. Singularities in 
Borel space are spaced by (exactly) |A|, see 0, while in the x plane their mirror 
singularities are, up to log corrections, as we see from regularly spaced by |A| _1 . 

Part of the results rely on a general property that if an analytic function is 
superexponentially small along a ray, then it becomes superexponentially large 
along a complementary ray. 

2. Main results 

In [2] it is shown that there is a unique analytic solution Fq of the associated 
equation 

m 

(6) ^^ (fe) =A(0,F ) 

k=l 

with ^ _1 Fo(^) — > 1 as £ — > 0. Fq has the property that 

(7) y(x) = F Q (x a e- x ) + 0(x- 1 ) 
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as x — ► oo if £ = x a e~ x stays bounded, and if £ s is a singular point of Fq, then y 
is singular along the array From JJJ) and the argument principle, if Fq(£) has 
a pole or an algebraic or logarithmic branch point at £ s then the actual solution y 
of has a quasiperiodic array of singularities, to leading order of the same type 
as that of Fq. The origin of the singularities of y(x) is the fact that in J5J the 
exponentials e~ kx become O(l) in the antistokes direction, and they can "pile up" 
to create blow-up of the solution. The singularities form 27rz-quasiperiodic arrays 
since they are governed by the size of x a e~ x . If y~ is the solution analytic in the 
lower half plane, then C~ = and S + can be calculated from J5J) and (0J . 
Assumptions. We recall A is analytic at (0,0). Let H — A(0, •) and $(s) = 
J S H(t)cLt. If to = 1 assume H is meromorphic and nonlinear. If to = 2 assume 
H is analytic except for isolated singularities and such that J a °° |$(s)| _1//2 ds < oo 
along a path P avoiding the zeros and singularities of <&. This includes the Painleve 
equations Pj and P/j, where H is a nonlinear polynomial 2 . Consider the one- 
parameter family T of solutions y of JQ) such that lim x ^ +oa y(x) = 0. 

Theorem 1. (i) Under the assumptions above, all functions in T except at most 
one are singular in the upper or lower half-plane. The singularities form arrays 
given by JfBjl . If H is polynomial of degree N > 1, then the local behavior of y near 
a singular point x s is to leading order [x — x s ) m ^ x ^ N \ 

(ii) If there exists anyo G T analytic in a sector of opening more than n centered 
on R + , then there is no other such yo, and S + = 0; yo is entire iff ya = A(x, 0) = 0. 

(Hi) For m = 1, if A(z, y) is polynomial of degree N > 1 in y and A(z, 0) ^ ; 
all solutions in T except at most one have infinitely many singularities in H^ +t , 
even if A(0, •) is linear. The singularity type is as in (i). 

2.1. General decay estimates. We state separately some useful and relatively 
basic estimates, helpful in proving Theorem ^ (hi). 1 

Lemma 2. If the function f ^ is analytic in a half plane {z : 'Si(z) > /?} and for 
all a > we have f(t) — 0(e~ at ) for t £ R + , t — » +oo, then for any e small 
enough f(z)e~ bz is unbounded in the closed right half plane for all b > 0. 

Lemma 3. Assume F S L 1 (IR + ) and for some e > we have 

(8) CF(x) = 0{e~ ex ) as x -> +oo 

Then F = a.e. on [0,e]. (The result is sharp as discussed after the proof.) 

Corollary 4. Assume F G L 1 and CF = 0(e~ ax ) as x — > +oo for all a > 0. Then 
F = a.e. on M+ . 

Corollary 5. Assume a > ^, / ^ is analytic in S = {z : \z\ > R, 2| argz| ^ 7r/a} 
and for any a > we have f(t) = o{eT at ) as t — > +oo. Then for any b 7 p > and 
< a' < a, the function f(z)e~ bz is unbounded in S . 

3. Proofs 

3.1. Proof of Theorem^ (i) Let first to = 1. We show by contradiction that Fq 
in © is not entire. Assume g(t) := Fo(e t ) was entire. We have g 1 = H(g). First, 

1 Lemma 121 resembles Carlson's Lemma 1111 pp. 185, but does not appear to imply it or to 
follow from it. Also, the referee points out that similar estimating methods appear in [§]• 
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H must have roots or else G = 1 /H would be entire and then the function Q given 

by 

(9) Q(g) = J G(g)dg + C = t 

would be entire, and injective since Q has an (entire) inverse. Thus Q is linear. 

Assume, without loss of generality, that H (0) = 0. Next we show that if H(t) = 
then g(to) ^ r for all to £ C. Indeed, if not, g'(to) — and we would have 

(10) fr<f)-ry = H(9U» -H(t) ^ 

v ' g(t) - t g{t) -r W 

But this is a contradiction, since to is a zero of finite order of g(t) — r and the left 
side of fTTijl tends to infinity. Since g(t) / for all t e C, we have g = e h with h 
entire, and since g is nontrivial it cannot avoid any other value, thus is the only 
root of H. Then ti = e- ,l H{e h ) where now the right hand side is meromorphic 
and everywhere nonzero, a case we have already analyzed. 

Let now m = 2, git) — Fo(e~*). We have g" = Hig). Then, multiplying by g' 
and integrating once, we get g' — + C. Using the condition g ~ e _t for large 
t and the analyticity of H for small g we get C = 0. Choosing a determination of 
the log and a to such that g (to) is not a zero or a singularity of $ we get 

/•s(t) 

(11) J(g(t ),g(t)) := / $(s)- 1/2 ds = t - 1 



's(to) 

We take the path P towards infinity. We have that J = lim^oo J{g{to), <?), .9 € -P, 
is finite and i s := fo + J is clearly a singular point of g. 

(ii) It follows from (@J that unless = every solution develops singularities 
in the upper or lower half plane (in both, except when C_ = or C + = 0). 
Conversely, if S + = 0, it follows from j2| that By is analytic in the C\K~ (Theorem 
1) and exponentially bounded in distributions in i?2u--e (Theorem 2). Thus CBy 
is a solution of the original differential equation with the asymptotic expansion y 
in any sector of opening less than 2ir + ir/2 centered on R + . In particular, this 
solution cannot be entire unless it is identically zero, and then Aq = 0. The type 
of singularity follows from and the discussion at the beginning of 33 

For part (iii) we need the following result. 

Lemma 6. Assume y solves 

N 

(12) y' = ^S fc (x)y fe 

k=0 

N > 1, in a sector S where the analytic coefficients Bk are O(o; 9fe ) and Bm ^ 0. 
Assume furthermore that for a sequence Xj — > oo in S, yixj) grows faster than 
polynomially in Xj. Then for all j large enough, within a distance Oix~ qN yixj) x ~ N ) 
of Xj there is a singularity of y. 

Note. Heuristically, near a large Xj, the dominant balance is of the form y' ~ 
Cy N x qN , 0, which forms a singularity as described. 

Proof. Let x — Xj + a and r\ = 1/y. The equation for a is 
da /- 2 

(13) % = ~ ^ kl ( i r GM 

ar] Efc=o^ h B N -k{x 3 +a) 
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with the initial condition o~(l/y(xj)) = 0. It is straightforward to check that in 
our assumptions, for large enough Xj, G(a,n) is analytic in the polydisk {(a,n) : 
H < 1,\t}\ < 2/\y{xj)\} where furthermore \G(a,rf)\ = 0(v N - 2 xj qN ). Then, 
in the disk {77 : \r)\ < 2/\y(xj)\} there exists a unique analytic solution a and 
cr = 0(w N ~ 1 xJ qN ). It is clear that Xj + er(0) is a zero of r\ and thus a singularity 
ofy. □ 

Proof of Theorem^ (Hi) Assume that there is a solution 7/0 of such that yo — > as 
x — ► +00 with finitely many singularities in a sector S of opening 7r + 2e centered on 
R + . (There may be an exceptional solution, for instance when the Stokes constant 
is zero and the asymptotic series converges.) Let yi be another solution of y' = 
A(x~ l , y) such that 7/2 —> as x — > 00 and let 5 = 7/2 — Ml- The equation for 5 is of 
the form 112(1 with Bq = 0. From the general theory of differential equations ^2] 
or from 3 it follows that S = 0(e~ x ) as x — > 00. By LemmaEJthere is a sequence 
of xj — > 00 in 5* so that S > const. exp(|a;j| 1 ~ £ ). Near every Xj with j large enough 
there is, by Lemma El a singularity of S and thus of 7/2- 

3.2. Connection with integrability. It is seen that if H is a polynomial of degree 
iV>2ifm=lorJV>3ifm = 2a one parameter family of solutions forms 
arrays of branch-point singularities in C and the Painleve test fails. For other 
values of p, the corrections Fj,j > 1 |2] can be calculated to determine the exact 
type of singularity. More can be done however. Assuming we are in the case where 
singularities are branch points, in view of the exact description of type and location 
of singularities it is possible to calculate the monodromy of solutions along a curve 
in C winding among sufficiently many singularities to show dense branching, a 
concept proposed by Kruskal, which can be used to show absence of continuous 
first integrals as in 0, |o]. This will be the subject of a different paper. 

3.3. Proof of Lemma l3l We write 

(14) / e~ px F{p)dp= / e~ px F(p)dp+ / e - px F(p)dp 

JO JO Je 

we note that 

e- px F(p)dp ^e~ eX J t \F(p)\dp ^ e'^WF^ ^ 0(e-' x ) 

Therefore 

(16) g{x) = [ e~ px F(p)dp = 0(e~ ex ) as x -> +00 

Jo 

The function g is manifestly entire. Let h(x) = e ex g(x). Then by assumption h is 
entire and uniformly bounded for i£l (since by assumption, for some xq and all 
x > xq we have \h\ ^ C and by continuity max \h\ < 00 on [0, xq]). The function is 
bounded by H-FHi for x £ iM.. By Phragmen-Lindelof 's theorem (first applied in the 
first quadrant and then in the fourth quadrant, with /3 = 1, a = 2) h is bounded in 
the closed right half plane. Now, for x = — s < we have 

(17) e— f\ sp F(p)dp^ f\F{p)\ < HFlIx 

Jo Jo 

Again by Phragmen-Lindelof (and again applied twice) h is bounded in the closed 
left half plane thus bounded in C, and it is therefore a constant. But, by the 
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Riemann-Lebesgue lemma, h — > for x — is when s — > +00. Thus h = 0. There- 
fore, with the characteristic function of A, 

(18) f F{ P )e- i ^dp = fi{ X[Qt€] F)=Q 
Jo 

for all s £ R entailing the conclusion. 

Note. In the opposite direction, by Laplace's method it is easy to check that for 

any small e > we have Ce~ p ' — o (e x J an d for any n C ^e~ En+1< - 1 ^ p ^ — 

o{e~ x / Lrl ^) where E n is the n-th composition of the exponential with itself and 
L n is the n-th composition of the log with itself. □. 

3.4. Proof of Proposition [2J By a change of variable we may assume that (3 = 
— 1. Assume that for some b > e~ bz f(z) was bounded in the closed right half 
plane. Then -ip(z) — (1 + z)~ 2 e~ hz f(z) satisfies the assumptions of Lemma |S] 
But then ip( z ) = CC~ 1 t\}{z) satisfies the assumptions of Corollary 01 and ip = 0. 

Note 7. There is indeed loss of exponential rate: the entire function T(x)^ 1 = 
0(e~ ninn ) on R + but is bounded by 0(e' a: ' 7r / 2 ) in the closed right half plane. 

4. Appendix: overview of Laplace transform properties 

For convenience we provide some standard results on Laplace transforms. 

Remark 1 (Uniqueness). Assume F £ L 1 (R + ) and CF = for a set of x with an 
accumulation point. Then F = a.e. 

Proof. By analyticity, CF — in the open right half plane and by continuity, for 
s £ I, CF(is) = = TF where TF is the Fourier transform of F (extended by 
zero for negative values of p). Since F £ L 1 and = TF £ L 1 , by the known 
Fourier inversion formula 0, F = 0. □ 

Lemma 8. Assume that c and f(z) is analytic in H c := {z : !ft z c}. Assume 
further that g(t) := sup c ,^ c |/(c' + it)\ £ L X (R, dt). Let 

(19) F(p) = — \ e^f(x)dx =: (C' 1 f)(p) 

/>OG 

Then for any x £ {z : 5ft z > c} we have CF = / e~ px F(p)dp = f(x) 

Jo 



Note that for any x' = x[ + iy[ £ {z : 5R z > c} 



(20) / dp 



C+lOO 



c— 'too 



e p(-x ') /(s) d | s |^ / ,00 dpe P(c-x' 1 )|| fl || 1< _M 1 



and thus, by Fubini we can interchange the orders of integration: 



1 



(21) U(x') = / e~ px — / e px f(x)dx 



l) 27TI J c _ ioo 



c+ioo 



1 rc+ioo i-OC 1 rc+ica f( T \ 

1 dxf(x) / dpe' px ' +px = — ; / -^-dx 
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Since g € L 1 there exist subsequences {r„},{— t^} tending to infinity such that 
\g{T n )\ — > 0. Let x' > 5Rx = xi and consider the box B n = {z : $lz <E [xi, a;'], £ 
[— r^,r„]} with positive orientation. 

(22) f l^-ds = -f(x') 

while, by construction, 

( 2 3, , im / iii* . m-i. - r~ ^ 

n ^ooJ Bn X'-S J x '- io o X'-S J c -ioo X'-S 

On the other hand, by dominated convergence, we have 
(24) / -y^-ds ^0 as x' -» oo 

J x' — ioo X — S 
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